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Abstract 

We use a perturbatively derived effective field theory and three-dimensional lattice 
simulations to determine the longest static correlation lengths in the deconfined QCD 
plasma phase at high temperatures (T>2T C ) and finite densities (/x<4T). For van- 
ishing chemical potential, we refine a previous determination of the Debye screening 
length, and determine the dependence of different correlation lengths on the number 
of massless flavours as well as on the number of colours. For non- vanishing but small 
chemical potential, the existence of Debye screening allows us to carry out simulations 
corresponding to the full QCD with two (or three) massless dynamical flavours, in 
spite of a complex action. We investigate how the correlation lengths in the different 
quantum number channels change as the chemical potential is switched on. 
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1 Introduction 



With the advent of RHIC and ALICE, there is a growing need for a precise understand- 
ing of various properties of QCD at temperatures up to a GeV. At the moment, we 
are still far from a satisfactory level in this respect, even for the equilibrium properties 
of the plasma. Indeed, even though the system can in principle be described as a gas 
of quarks and gluons, a fully perturbative computation with these degrees of freedom 
does in practice not work well at any reasonable temperatures below ~ 10 10 GeV, since 
the perturbative series is badly convergent due to infrared sensitive contributions fl 
[|J. On the other hand, the only systematic fully non-perturbative method available, 
four- dimensional (4d) lattice simulations, is severely restricted in the presence of light 
dynamical fermions, and even more so in the presence of a finite baryon density; the 
current state of the art is summarised in || . A first principles solution of any problem 
related to non-equilibrium phenomena remains even further in the future. 

The physics problem we address in this paper is the determination of various static 
correlation lengths. To tackle this problem, we will employ a method which can over- 
come some of the difficulties mentioned above. The method combines analytic and 
numerical techniques, in a way that both are used in a regime where they are well 
manageable. First, perturbation theory is employed in deriving, by what is called di- 
mensional reduction || [7|, ||, || |3], £|] , an effective theory for the long range degrees of 
freedom of the system. At high temperatures all such degrees of freedom are bosonic, 
since fermions are screened by non-zero Matsubara frequencies irTn, with n odd. The 
reduction step can thus be carried out with dynamical, massless fermions and, as we 
shall see, with a finite chemical potential. Second, non-perturbative lattice simulations 
are used to study the infrared sensitive dynamics of the remaining degrees of freedom. 
The theory to be studied with simulations is the SU(3)+adjoint Higgs model in three 
dimensions (3d), and many of its properties have already been determined [iTcfl . 

The method we employ suffers clearly from a number of restrictions. First of all, in 
QCD it is limited to temperatures above T c , roughly T > 2T C []. This can be inferred 
from a comparison of the dynamical scales described by the effective theory to those 
integrated out [|J] as well as, more concretely, from a direct quantitative comparison of 
the correlation lengths measured within 3d for SU(2) ]T3| , with those determined using 
4d SU(2) lattice simulations at finite temperature (without fermions) |14]Q- Second, as 
we shall see, it is also limited by the largest chemical potential that can be reached. 
We will be able to go up to /i < 4T. 

If there are restrictions, there are also strengths. What can be done within the 
effective theory can be done quite precisely. The simulations are comparatively easy 

1 For the electroweak theory, in contrast, similar methods allow for a precise determination of the 
properties of the phase transition at T = T c for Higgs masses m# ~ 30. ..250 GeV |ll], |i~2| . 

2 For early work with similar conclusions both in SU(2) and SU(3), see Q. Recently, analogous 
results have been reached also by considering gauge fixed correlators |L{|, as well as by considering 
dimensional reduction of pure SU(3) from (2+l)d to 2d Jl6| . 
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technically, they yield high accuracy continuum limits for the quantities under investi- 
gation, and hence there is little ambiguity in the conclusions. For instance, studies in 
3d have produced detailed and accurate insight into the relative sizes of perturbative 
and non-perturbative contributions to the Debye screening length [L3|], a quantity 
which could be directly relevant for such signals of the quark-gluon plasma as J/\I/ 
suppression. 

In the present work, we extend the study of 0] first to SU(3) Yang-Mills theory, 
and then to two- and three-flavour QCD in the massless limit. In comparing with 



earlier results in SU(3) 0, [17[], we improve on the accuracy of Debye screening length 
determination and measure the screening lengths also in the other quantum number 
channels. We also discuss explicitly the N c and Nf dependence of our results — here 
N c is the number of colours and Nf is the number of (massless) quark flavours. 

We then proceed to extend our calculations to finite baryon density. Our approach 
is straightforward: dimensional reduction leads to an additional complex term in the 
effective action, which has to be absorbed in the observables for practical simulations. 
This "reweighting" displays in principle similar problems as chemical potential simula- 
tions in four dimensions (for reviews, see |5|, M). We find that for the effective theory 
there exists a range of volumes and ratios /x/T, however, for which the problems are 
in practice manageable. We demonstrate this by investigating a range of imaginary 



chemical potentials [19], for which we find complete agreement between reweighted cal- 
culations and those using the exact action. We then employ the reweighting technique 
to the case of real chemical potentials. 

Let us note that our approach relies on ix/{txT) being small, but no restriction is 
imposed on fermionic masses m;, i — 1, Nf, which we take to be zero. In 4d simula- 
tions, on the other hand, progress has been possible when the m; are taken to be very 
large E0, |2"T], whereas the temperature can be small. 



2 Continuum formulation for \i = 

We start by reviewing the result of the dimensional reduction step for fi = 0. We follow 
closely Ref. @. 

The effective theory emerging from hot QCD by dimensional reduction || is the 
SU(3)+adjoint Higgs model with the action 

S = | C / 3 x{^Tr4 + Tr[A,Ao] 2 + m2Tr^ + A 3 (Tr^) 2 |, (2.1) 

where F^ = diAj — djAi + ig[Ai, Aj], Di = di + ig^Ai, Fij, Aj, and A are all traceless 
3x3 Hermitian matrices (Aq = A^T a , etc), and g\ and A3 are the gauge and scalar 
coupling constants with mass dimension one, respectively. The physical properties of 
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the effective theory are determined by the two dimensionless ratios 

I= * y= ^ 3 = 9i) : (22) 

where /2 3 is the MS dimensional regularization scale in 3d. The parameters well 
as the scale g 2 , are via dimensional reduction functions of the temperature T and the 
QCD scale A^. For N c = 3, the result is || 

2 _ 247r2 1 

93 " 33 - 2N f Hf, g /A m ) ' (2 ' 3) 

x = l^Ii I (24 ) 

33-2JV) Hfi x /A m Y { • ) 

(9-Nf)(6 + N f ) 486 - 33N f - UN 2 - 2N 3 f 

»<*> = 144A + 96(9 -A^ i + OW ' (2 ' 5) 

We have assumed here Nf massless flavours, but the inclusion of quark masses is also 
possible in principle. 

The scales entering in Eqs. (|2.3|), (|2.4|) are 



\ii = iirfah (2-6) 
fi'T = 47re-»T « 7.0555T, & = exp (Z^t^^l^ ( 2 .7) 

54- 22iV, 4 . . 

c 9 = 1, c x = 9 _ N / + - 3 N f . (2.8) 

To be explicit, for Nf = 2 we obtain 

2 _ 24?r 2 T 7 _ 7 15 

93 ~ 29^(8.11^^)' X ~ 291n(6.91T/A 5J g)' V ~ ISA + 28^' ( ^ 

Corrections to these expressions are of relative magnitude ~ 0(a s /ir) 2 ~ 0(x 2 ). 

An important question is how small one can in practice take T/A^jg and still have 
only small corrections in Eqs. ( |2.3| )-(2T5). As was observed in for y(x) the expansion 



appears to converge surprisingly fast, with an error on a few percent level even at 
T ~ Ajyfg. The expansions need a priori not be as good for every parameter, however. 

Consider for instance g\ . One may expect higher loop graphs to amount to an in- 
crease in the effective scale factor j2 g , since there are many "massive" (~ 2ttT) particles 
inside the loops instead of one. Thus an upper bound on the error in g\ can be ob- 
tained by evaluating the 2-loop QCD running coupling at the 1-loop scale £L g ~ 7T 
(for Nf = 0). This gives a correction of up to 25%, which is quite large. An explicit 
computation of g\ including effects of order 0(g 6 ) would thus be welcome in order to 
clarify whether the expansion for g 2 in reality converges as rapidly as for y(x) or not. 



3 



Such computations are beyond our scope here and we will simply use Eqs. (|2.3|) - (|2.5|) ; 
a comparison with direct 4d simulations at Nf = turns out to be quite satisfactory 
within this procedure, and thus we will work under the assumption that the error in 
<7§, x is similarly small as in y. 

Finally, we may sometimes want to express the temperature to which our simulations 
correspond in terms of T c , rather than A^jg. In view of Eq. ( |2.9| ), this requires knowledge 
of the deconfmement temperature T c /Aj^g, which can only be fixed with 4d simulations. 
For the case N f = 0, we take the value T c /A m = 1.03(19) from For N f °> the 

determination of T c is not accurate, and we will take T c /A^ = 1.0 as a reference. Our 
results however only depend on T/Aj^g, so they can be reinterpreted as corresponding 
to some other temperature relative to T c if need be. 



2.1 Operators and quantum numbers 

The physical observables on which we shall focus in this paper are spatial correlation 
lengths of QCD at finite temperatures. For the practical calculations we reinterpret 
the 3d theory in Eq. Q2.1D which we use to compute them, as a (2+l)d one. We 
take operators to live on the (x\, X2)-plane, and the correlations are taken in the x^ 
direction. Thus we compute the spectrum of the 2d Hamiltonian of the effective theory, 
whose eigenvalues then correspond to screening masses or inverse spatial correlation 
lengths of the 4d theory at finite temperature. Thus, unless otherwise stated, "bound 
states", "glueballs" etc. refer to the eigenstates of the effective Hamiltonian and to the 
screening masses of the 4d theory in the above sense. 
We use the quantum number notation 

(2.10) 

A, -A), (2.11) 

(2.12) 

The action in Eq. (|2.1| ) is invariant under these operations. In the finite temperature 
context, the operation R in the 3d theory is a remnant of the 4d time reversal oper- 



R : 


A - 


->D i} 


p. . _ 


* Fiji Aq — ► — Aq, 


P : 


A - 


-A, 


A - 


-* — A, F u — ► — F u 


C : 


A- 




p. . _ 


-> ~F*j, A — > -A* . 



ation ||24|| . The parity P means a parity on the 2d plane, i.e. a reflection across the 
Xi-axis. The charge conjugation C is a non-trivial quantum number only for SU(3), 
since for SU(2) it is just a global gauge transformation ir 2 , so that there are no gauge 
invariant operators odd in C |2"4| . In addition to these discrete transformations, we 
define a rotation in the (xi, X2)-plane, with the corresponding angular momentum J. 
Thus the full symmetry group is 5*0(2) £g> Z 2 (R) ® Z 2 (P) <8> Z 2 (C), and accordingly we 
classify our operators and states by . 

We note from Eqs. (|2.10|) — (|2.12|) that apart from A, the action of RPC on the 
operators in the (xi, X2)-plane corresponds to complex conjugation (and X2 — > — x 2 ). 

3 A recent computation p3| favours a slightly larger central value, but well within the error bars. 
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Thus any real operator which does not contain D2 is even in RPC . This means, in 
particular, that for the scalar J = operators in which no single D2 can appear, 
there are effectively only two quantum numbers left out of the original three, thus four 
different channels. The lowest dimensional gauge invariant operators of these types in 
the J = channels are: 



PC 
R 

PC 
R 

PC 
R 

PC 
R 



C++ 

or 

o: + 
o + - 



Tr^Tri^,... 
TrF? 2 ,Tr A 2 F 12 , 
TrA F 12 ,... 
TrA 3 ,TvA F? 2 ,. 



(2.13) 



All operators here with C — — 1 vanish identically for SU(2). Going out of the plane 
(i.e., allowing also for F 13 , F23), other channels would in principle become possible |24[ 
but we do not expect them to change our conclusions. 



3 Lattice formulation for fi = 

To simulate the effective theory in Eq. ( |2.1| ) on the lattice, we use the discretised version 
S 

X < i ■> -■> - < , 

X,i>j V 7 

- Tr (^(x)(/ i: (x) V (x + ^l(x)) + A^(2Tr (^(x)^(x)) - l) 2 , (3.1) 



l-TTReTr^(x) +2^Tr(^(xMx)) 

X,i>J J x 



where the fields in the lattice action have been rescaled relative to the continuum, 
v4o(x) = (ft/a) 1//2 <£>(x) where a is the lattice spacing, and ^ij(x) denotes the elementary 
lxl plaquette in the i, j-plane located at x. The parameters of the continuum and 
lattice theory are up to two loops related by a set of equations which become exact in 
the continuum limit [p5 , 26|] , 
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+ 



2 ' 


A = 








- -3 


18 \K 


1 

16tt 2 


(60x 



3xk 2 




3.1759115/5 
4tt 



(3.2) 



For a given pair of continuum parameters x, y these equations determine the lattice 
parameters n, A as a function of the lattice spacing, and hence govern the approach to 
the continuum limit, (3 — > oo. We have not implemented 0(a) improvement [27j here, 
since the discretisation effects in the correlation lengths are already quite small at the 
values of (3 we are using. 
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3.1 Observables 



Operators for all the quantum number channels considered can be constructed from a 
number of basic operator types. We consider operators involving only scalar fields or 
products of scalar and gauge field variables, 

i? 2 (x) = Tr(^ 2 (x)), 

i? 3 (x) = Tr(^ 3 (x)), 

Li(x) = TV (^(x^xMx + i)U} (x)), 

5x(x) = Tr(y?(x)f%(x)), 

5 2 (x) = Tr(^ 2 (x)^(x)). (3.3) 
Furthermore, we have loop operators constructed from link variables only, 

Cj xl (x) = Tr (c/ i (x)C/ i (x + ?)C/t(x + i)C/j(x)) , i,j = 1,2, i ^ j, (3.4) 

and in addition to the elementary plaquette C lxl , we also consider squares of size 2x2 
as well as rectangles of size 1x2, 1x3, 2x3. Another useful pure gauge operator is 
the Polyakov loop along a spatial direction j, 

L-l 

Pf >(x) = ReTr J] ^(x + mj), j = 1, 2, (3.5) 

m=0 

which can be used to extract the 3d string tension. It also gives useful information 
about finite volume effects via so called torelon states [^, ^J. 

Operators with definite quantum number assignments are constructed from the 
above types by taking linear combinations with appropriate transformation proper- 
ties. Clearly operators containing an even number of scalar fields, such as B 2 , will 
couple to R — +1 states, and those with odd, such as Bi, to R — — 1. The different 
P, C channels can be chosen by utilising the projection operators \{1 ± P), |(1 ± C) 
where, e.g., C-Bi(xi,x 2 ) = — Bl(xi,x 2 ), PBi(xi,x 2 ) = B*(x\, — x 2 ). Different spin 
states are obtained by employing the operation R{0 n ) which rotates the operators by 
n(7r/2) around x. Note also that even though lattice rotations are restricted to mul- 
tiples of 7r/2, we are in practice close enough to the continuum limit that continuum 
symmetries are reproduced within the statistical errors; thus we keep the continuum 
notation in terms of J even on the lattice. 

In the following we list the operators used corresponding to Eq. ( |2.13|) : 

channel: 
R: i? 2 (x), 

L: Re(L 1 (x) + L 2 (x)), 

C: real part of symmetric combinations of C lxl , C 2x2 , C lx2 , C lx3 , C 2x3 , 
P: P 1 (i) (x) + P 2 (i) (x), 
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T: (Pf'(x)) ! + (#(x)) ! , 
channel: 

C: imaginary part of symmetric combinations of C lxl , C 2x2 , C lx2 , C lx3 , C 2 
B: ImEn=i^n)5 2 (x), 

0l + channel: 

£: ImEn=i^n)£i(x), 

0l~ channel: 
i?: i? 3 (x), 

5: ReE^Uii2(0n)5i(x). 



We have also measured higher spin states. For 2+ + we have a fairly large basis, identical 
to the one described in |13|]. The J = 1 states prove to be quite heavy, so that their 



relevance for the 4d finite temperature system is not immediately obvious at modest 
temperatures, and thus for simplicity we consider here only the ground states in the 
channels with R = ±1, without specifying the other quantum numbers. 

3.2 Blocking and matrix correlators 

For later reference, let us recall the general principles of how the operators discussed 
above can be used for a reliable extraction of correlation lengths . 



The eigenstates of the 2d Hamiltonian in the region of parameter space where we 
work are bound states. In order to increase the overlap of our operators onto such 
extended states, we construct smeared link and scalar field variables fa as described 
I3fl . For every quantum number channel, we then measure the correlation matrix 



in 



C l3 {t) = (0l(t)0,(O)), (3.6) 

where we have denoted t = x 3 . This matrix can be diagonalised by solving a generalised 
eigenvalue problem, 

C7- 1 (0)C(t)v„, = A n (t)v n , (3.7) 

where \ n (t) ~ exp(— aM n t). We carry out this procedure at t = a, 2a, 3a, which gives 
somewhat differing eigenvectors v n , and check that the final outcome remains the same 
within error bars. We normalise the eigenvectors according to vj l C(0)v ri = 1, so that 
= Ei v nVi satisfy $^$ n = 1, and are thus the normalized eigenstates of the 2d 
Hamiltonian. The final results are extracted from the correlation functions 

G„(t) = (*i(t)*„(0)) = vtC(t)v„, (3.8) 
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by computing effective masses 



aM eSjn (t) = - In 



G n {t+1) 
G n (t) 



(3.9) 



and ensuring that these have attained a stable plateau value. Information about the 
composition of $ n in terms of the operators used in the simulation, which we shall 
quote in some of our tables, is obtained from the overlaps (0j$ n ) = J2j Cy(0)v^. 



3.3 Simulation and analysis 

In our Monte Carlo simulation of the lattice action, Eq. (|3.1|) , link variables are updated 
by a combination of heat bath and over-relaxation steps with algorithms described in 
|3T |. The scalar fields are generated by a combination of heat bath and reflection 
steps [ 3"2"| . One "compound" sweep consists of several over- relaxation and reflection 
updates following each heat bath update of gauge and scalar fields. Measurements are 
taken after every compound sweep. Typically, we gathered between 5 000 and 20 000 
measurements depending on the lattice sizes. Statistical errors are estimated using a 
jackknife procedure with bin sizes of 100 - 250 measurements. 



4 Numerical results for \i = 

In this section we discuss simulations of the effective theory for parameter values cor- 
responding to hot SU(3) gauge theory with Nf = 0, 2, 3, 4 flavours of fermions at zero 
baryon density. Detailed investigations of finite volume effects and scaling of correlation 
lengths have been performed for the 3d pure SU(2) and SU(2)+adjoint Higgs theories 
in [L3J, as well as for the 3d pure SU(3) theory in |31], [17| . In these works explicit 
continuum extrapolations were performed. For large enough (3 the scaling behaviour 
was found to be quite good, with continuum results differing from those on a fine lat- 
tice by only a few percent. Our gauge couplings here were chosen to produce a similar 
lattice spacing, and checks also show good scaling. The parameter combinations and 
the lattices used for them are collected in Table [p. 

4.1 The mass spectrum for Nj = 

The main features of the spectrum, displayed in detail in Tables |5], |7| in the Appendix, 
are the same as those found in previous studies of the confinement region for SU(2) 



with fundamental [29] or adjoint [13] scalar fields. There is a dense spectrum of bound 
states, consisting of a replication of the glueball spectrum found in the d = 2 + 1 pure 
SU(3) theory and additional bound states of scalars, with little mixing between 
the two. This may be concluded from the fact that, as indicated in Table |5], some 
eigenstates are composed predominantly of purely gluonic operators C, with practically 
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0.09191 


0.4830 


L 2 ■ 


T = 


30 3 


L 2 ■ 


T = 40 3 
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1.5 


0.08680 


0.47890 


L 2 ■ 


T = 


30 3 
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0.07814 


0.52741 


L 2 ■ 


T = 


30 3 






4 


2.0 


0.062921 


0.569680 


L 2 ■ 


T = 


30 3 







Table 1: The lattice parameters and sizes used for calculations at // = 0. Let us stress 
that for Nf > 0, T/T c means really T/Ay^, see Sec. [|. 



tPC 
J R 


gauge-Higgs 


pure gauge 


scalar glueball 


glueball 




0.994 (19) 2.511 (65) 


2.575 (18) 




2.95 (16) 3.61 (32) 


3.841 (28) 




2.46 (10) 3.50 (25) 


3.795 (27) 


2| + 


3.355 (94) 4.14 (28) 


4.257 (33) 



Table 2: Mass estimates M/g 2 atT = 2T C and Nf = 0. Our data for the glueballs are 
compared with the results obtained in the pure gauge theory 



no contribution from operators carrying the same quantum numbers but containing 
scalars. Comparing with ||3l| , we find that our gluonic states have quantitatively 



the same masses as the corresponding glueballs in the pure Yang-Mills theory; the 
comparison is shown in Table |^. 

As can be observed from Table [|, however, in the physically interesting region of 
couplings the lightest states in given quantum number channels are not glueballs but 
always states including the scalar field tp. This can be understood in the sense that 
in the finite temperature context the phase transition is (for Nf = 0) assumed to be 
driven by the Z(N) symmetry related to A (or if on the lattice), and thus if is the 
dominant infrared degree of freedom. Let us now discuss the corresponding correlation 
lengths in more detail. 



4.1.1 The "magnetic" sector, R = +1 

Consider first operators with R = +1. We will call this the "magnetic" sector, even 
though the operators can include an even number of Ao-fields. This sector determines 
the correlation length related to the lowest lying glueballs and all other 0^ + observables, 
as well as that felt by the real part of the 4d temporal Polyakov loop |£| Q (not to 
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Figure 1: The spectrum of screening masses in various quantum number channels at 
N f = 0,T = 2T C (left), N f = 0,T ~ 10 n T c fhp/it/ Fz/fed sym&o/s denote 5c? ^«e6aiZ 
states, which have become the lightest excitations at T ~ 10 n T c . T/ie states 1 + are 
much heavier than 1_ at high temperatures, and thus not visible on the right. 



be confused with the spatial Polyakov loop in the 3d theory). In order to express 
the masses in units of temperature, we use the perturbative expression for the gauge 
coupling in Eq. Q to arrive at gl(T = 2T C ) « 2.7T and g 2 (T = 10 n T c ) « 0.25T. 
The spectrum in these units is shown in Fig. |l|. 

Let us first discuss the channel 0^ + . From Table |5| in the Appendix, we note that 
the lightest state (open circle in Fig. |l|(left)) is dominantly ~ TrAg, while the next 
state (filled circle in Fig. |l](left) ) is dominantly ~ Tr F^ 2 . Thus, the naive ordering of 
states consisting of the fields Aq ~ gT and A, ~ g 2 T is reversed non-perturbatively: 
at low temperatures it is the Ao which is responsible for the lightest gauge invariant 
state in the system. 

As the temperature is increased, however, the ordering gets changed. The lowest 
state becomes heavier with increasing temperature, as the scalar "constituent mass 
parameter" y 1 ! 2 is then growing (for a discussion of the constituent picture in a similar 
context, see |34] , |15|]). On the other hand, the glueball ~ TrFf 2 does not contain 



much Aq admixture. Correspondingly its mass, in units of o|, is quite insensitive to y, 
so that for large enough temperatures indeed the glueball corresponds to the lightest 
excitation. This situation is, however, only realised at high temperatures T > 10 2 T C ; 
an extreme example is shown in Fig. pjright). In this limit the Aq may be integrated 
out, leaving the 3d pure SU(3) gauge theory as an effective theory ||. 

For the channel 0T~, the behaviour is similar, but the masses are much larger. Let 
us mention here already that in the 4d theory with fermions at modest temperatures 
T ~ 2T C , all the higher-lying bosonic states in general are probably heavier than some 
gauge invariant states consisting of fermionic fields ~ ijjijj, which are not addressed by 
our theory, and whose mass we may expect to be < 2nT. (At lower temperatures such 
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states can be even lighter, since in the chiral limit they are expected to represent the 
critical degrees of freedom.) 

To get an impression of the quality of dimensional reduction for Nf = 0, we now 
compare the 0+ + state ~ Tr^Q with that found directly from a 4d simulation at 
Nf — |4j. Our value is M[0++] = 2.71 (6)T, to be compared with M = 2.60(4)T 



reported in [jT3[| . In complete analogy with the SU(2) case |13|], it is thus found that 
dimensional reduction quantitatively reproduces the lowest screening mass of the 4d 
Yang-Mills theory at temperatures as low as T = 2T C . 

We can also carry out another comparison. The real part of the 4d temporal Polyakov 
loop carries quantum numbers 0+ + , and non-perturbatively mixes with other operators 
in that channel. Correspondingly, its decay should also be determined by the bound 
state Tr Aq. Indeed, the exponential decay of the 4d temporal Polyakov loop correlator 
is in 4d measured to be ~ 2.5T which is quite compatible with the value ~ 2.7T 
we find here. 

On the other hand, in contrast to SU(2) where even several excitations agree between 
the full and the effective theory, we find here deviations of about 20% in comparing 



with [|14j]. This may signal that the effective theory is not as reliable for states whose 
mass is approaching ~ 2nT. However, we also believe that there could be some room for 
improvement upon this apparent disagreement. The 4d simulations are quite difficult, 
and it is not easy to guarantee at this stage that the infinite volume and continuum 
limits are reached for the excited states. 



4.1.2 The "electric" Debye sector, R = — 1 

Of particular phenomenological relevance for the QCD plasma is the Debye mass, 
whose inverse gives the length scale over which the colour-electric field is screened. 
The expansion in powers of coupling constants for this quantity reads 



9l 



m 



LO 
D 



9l 



+ -^ln 

4:71 



jf 0(g 3 T) 
+ c Nc + 



The perturbative leading order result is |36 



m 



to 

D 



Nc 
3 



+ 



N i 



1/2 



9l 



gT = gl^+0(g 3 T), 



(4.1) 



(4.2) 



but at next-to-leading order ~ g%, only a logarithm can be extracted |37J (the 2nd 



term on the right hand side of Eq. fl4.1|) ), whereas the coefficient cn c is entirely non- 
perturbative. To determine it, we use the gauge invariant definition of the Debye mass 



based on the Euclidean time reflection symmetry as given in ||24| . According to this 
definition, the Debye mass corresponds to the mass of the lightest state odd under this 
transformation. The remnant of Euclidean time reflection symmetry in our reduced 
model is the scalar reflection symmetry R, and the Debye mass thus corresponds to 
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m D /gl 


pert, part 


c 3 


0(9 3 T)/gl 


m D /T 


T = 2T C 
T ~ 10 n T c 


1.70(5) 
3.82(12) 


0.514 
2.165 


1.65(6) 
1.65(6) 


-0.46(6) 
0.00(12) 


4.6(2) 
0.96(3) 



Table 3: The different contributions to the Debye mass, Eq. t \4-A )- The perturbative 
part means the 1st and 2nd terms on the right-hand- side of Eq. U-A )- 



the mass of the lightest R — — 1 state of our spectrum. As shown in Fig. [I], this is the 
0Z + ground state ~ Try4o-Fi2, while 0l~ ~ Tr is slightly heavier. 

With this definition, the coefficient cn c can be measured separately from the expo- 
nential decay of a Wilson line in a 3d pure gauge theory [E4|. The measurements for 



N c = 2, 3 have been performed in jl7| with the results = 1.14(4), C3 = 1.65(6). On 
the other hand, the mass we have measured here also includes the 0(g 3 T) correction 
in Eq. (|4.1|) . We can thus now estimate the magnitude of the remainder; our results 
are shown in Table [| We find that the 0(g 3 T) corrections are less than 30% even at 
temperatures as low as T = 2T C , and they disappear entirely for asymptotically large 
temperatures. The sign is negative, so that the complete result is somewhat smaller 



than estimated in [17| based on the perturbative contributions and C3 alone, ~ 5T 
instead of ~ 6T. The result is also about 20% smaller than what can be extracted 
from ||, a difference which we presume to be due to the relatively small lattice sizes 
used there. 

We can conclude that for temperatures of physical interest, the Debye mass according 
to this definition is entirely non-perturbative, with its next-to-leading order correction 



being larger than the leading term. This remains to be the case up to T ~ 10 7 T C [|17 
But even at T ~ 10 n T c , as the table shows, there are sizeable corrections to the leading 
behaviour. This leads to the conclusion that the scale dominating the Debye mass is 
~ g 2 T with non-perturbative physics for all temperatures of interest, in contrast to the 
naive expectation ~ gT. The reason for this behaviour is the large non-perturbative 
coefficient cn c that appears in front of the correction term in the series in g. 

4.1.3 iV c scaling between SU(2) and SU(3) 

Let us finally discuss the scaling with N c . In |yj it was found for SU(iV c ) pure gauge 
theories with N c — 2, ...5 that the differences in the mass spectra can be accounted for 
by the leading order 1/AT? corrections in a large N c expansion, and that the coefficients 
of these are remarkably small. In theories with various scalar fields, the glueball content 
has been found to be practically identical to that of pure gauge theories for N c = 
2, 3 PPI, BR O], and thus the scaling behaviour with iV c is preserved. It is then natural 



to ask if the same scaling behaviour holds for the scalar bound states, which are not 
present in the Yang-Mills theory. In Fig. |^ we plot some of our low lying states from 
the current analysis as well as the SU(2) case |L3] (we have to focus on states with 



12 



T=2T , N,=0 





_L 

-f- 






~m - 


- 9 




-e 


-o 

0'. states 



0.1 



1/n; 



0.2 



0.3 



1.4 

i 1-0 

z° 

S 1-2 
0.8 
0.4 



T=2T C , N f =0 



0.1 



- - ■ 

- - -□ 



□ 2 + states 
0~ + states 



0.2 



0.3 



1/n; 



Figure 2: The scaling with N c of some of the low lying 0+ + (left) and 0_ + , 2\ + (right) 
states. Filled symbols denote glueballs, and the lines are to guide the eye only. 



C = +1 in order to have an SU(2) counterpart). Indeed we observe a similar scaling 
for the scalar states as for the glueballs. This suggests that the screening masses of 
hot SU(3) gauge theory are close to the N c — > oo limit, and hence large N c methods 
may be useful approximations to analytically deal with some of the non-perturbative 
aspects discussed here. 



4.2 The mass spectrum for Nf ^ 

Having convinced ourselves that dimensional reduction of pure Yang- Mills theory works 
at the very least for the lowest static bosonic correlation functions of the theory, we 
may now move ahead to exploit the advantage of this formalism, the easy inclusion of 
fermions. For the case of massless fermion, all that is required is to evaluate Eqs. ( |2.3|) - 
( |2.5| ) for the desired Nf and simulate the effective theory with the corresponding pa- 
rameters. Each of these simulations gives results completely analogous to those of the 
Nf = case, except for a slight shift in the values. For the reader interested in the 
detailed numbers, we collect our data in Tables || |8| in the Appendix, at tempera- 
tures T = 2Ajyfg, and 1.5Aj^g, which probes the lower limit of the range of validity of 
dimensional reduction. 

With the tables at hand, the spectrum shown in Fig. [I] for Nf = is now also known 
for Nf = 2,3,4. In Fig. |^ we display the dependence on Nf of the low lying J = 
states. The general behaviour of slightly rising mass values with Nf can be understood 
from Eqs. ( |4.1| ), ( |4.2j ). Increasing Nf increases the value of m^,° corresponding to the 
bare scalar mass, and hence the masses of the bound states increase as well. In units 
of the temperature, the increase from Nf = to the phenomenologically interesting 
case Nf = 2. ..3 is about 30% for 0| + , < 20% for the other channels. The correlation 
lengths decrease accordingly. 
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Figure 3: The Nf dependence of the J = spectrum at T = 2T C (to be more precise, 
for Nf > T c means really A^g, see Sec. \^). Filled symbols denote glueball states. We 
should stress that only the part M < 2nT of the spectrum can be expected to directly 
represent the lowest states in 4d finite temperature QCD. 



4.3 The spatial string tension 

Finally, let us discuss one observable of the 3d theory which does not have an inter- 
pretation as a physical 4d correlation length. A Polyakov loop in a spatial direction 
couples to a flux loop state (or torelon) that winds around the periodic boundaries of 
the finite volume. The exponential fall-off of its correlator is related to the mass of the 
flux loop, according to 

]T (p} L \x + t3)Pj L \x)) ~ e- aMp{ - L)t , aM P {L) =a 2 a L L. (4.3) 

x,i=l,2 

Such a flux loop state can be easily identified through its energy scaling linearly with 
the size of the lattice as seen in, e.g., Table ||. Since the scalar fields are in the adjoint 
representation, they cannot screen the colour flux represented by the Polyakov loop 
in fundamental representation, in contrast to the case with fundamental scalars 



Thus, a string tension can be defined by the slope of the static potential at infinite 
separation, just as in pure gauge theory. Accordingly, the coefficient ol corresponds 
to the string tension at separation L. For large enough L, an estimate for the string 
tension in infinite volume is then provided by the relation [HO 



aVoo = a 2 a L + (4.4) 

We have extracted the string tension by diagonalising a basis of various smeared 
Polyakov loops, and the results are given in Tables |H]. 
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We note that the string tension in continuum units is rather insensitive to the precise 
parameter values y 5 clS it fluctuates by less than ~ 4% between the cases considered. 
Moreover, all values are similarly close to the ones obtained in pure gauge theory |H . 
This is a further indication of the insensitivity of the pure gauge sector to the presence of 
the adjoint scalar fields. In the context of Yang- Mills theory at finite temperature, the 
string tension of the effective theory corresponds to the spatial string tension as mea- 
sured in 4d simulations at finite temperature. Indeed, our value ^fo^jg 1 = 0.569(10) 
on the finest lattice for T = 2T C , Nf = is in good agreement with the spatial 4d result 
^W~ s /g 2 T = 0.586(5) given in gl 



5 Extension of the method to finite density 

As we have mentioned in the Introduction, there is no satisfactory algorithm to simulate 
lattice QCD in four dimensions at finite baryon density, with small quark masses. On 
the other hand, in the experimental situation of heavy ion collisions there always is a 
net baryon density, which in the collisions at and above AGS and SPS energies can be 
estimated to correspond to /x/T<4.0 [fE^] . As long as the temperature is sufficiently 
above T c , dimensional reduction is applicable for the lowest lying static correlation 
functions in QCD, as we have seen in the previous sections. We now proceed to 
dimensionally reduce QCD at finite temperature and density. As we shall see, it is 
possible to perform simulations in the /x/T-range of interest in this framework. 



5.1 The effective theory with /i ^ 



At leading order, the introduction of fi ^ leads to a very simple modification of 
the effective 3d theory in Eq. ( |2.1| ). New effects come only from fermions, where we 
change pf — > pf — ifi in the loop momenta; here pf denotes the fermionic Matsubara 
frequencies. The 1-loop effective potential computed in tells then that the effective 
action for the Matsubara zero mode of Aq has for fi ^ the terms (N c = 2, 3) 



V- 



A = 9 



~/iV c N f \ 
T 2 (— + ^-) + u 
V 3 6 > P 2tt 2 



N 



f 



Tr A 2 +ig 3 ti ±yfrTX - a3 1 " 4 



3tt 



24tt 2 



(TtAI) 2 . (5.1) 



Note furthermore 



There is of course no term linear in A , unlike in the Abelian case J44 
that for SU(2), Tr A 3 = 0. 

In general, the inclusion of /i 7^ also leads to other new operators than TrAg. 
However, as can be seen from the effective potential computed in [E3J, there are no 
higher order operators involving only A Q , at least at 1-loop and 2-loop levels. On 
the contrary, there could be operators such as ~ AqF^. As is usual in effective field 
theories, we expect such non-renormalizable operators to give contributions suppressed 
with respect to the dynamical effects arising within the effective theory in Eq. (|5.1| ) 
by a power of the scale hierarchy, and therefore we ignore them here. It is perhaps 
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Figure 4: The function T>(u) from Eq. (|5.4j ), compared with the small-u; limit. 

also worthwhile to mention that since no 75 appears in perturbation theory in QCD 
and since the construction of the effective theory in Eq. ( |5.1| ) is only sensitive to the 
ultraviolet and thus purely perturbative, no Chern-Simons type operators are expected 
to be generated. 

Going to 3d units (A* d T^A™, Jdrd 3 x T' 1 Jd 3 x) and denoting 



S z = Jd'xgiTrAi (5.2) 
the dominant changes due to /i 7^ in the action in Eq. fl2.1|) are then 



EEL. „^Ji + (JLV^El 



S^S + i zS„ z= ^ ; y ^ y[l + [ ^ )w:tWi y (5.3) 

Thus, one new operator is generated in the effective action, and one of the parameters 
which already existed, gets modified. We note that the new operator is quite special 
and it, for instance, does not generate new ultraviolet divergences for the parameters 
in the original action in Eq. ( |2.1| ). 

A more precise computation of the effective action requires the inclusion of 2-loop 
effects in Eq. (|5.1| ), as well as a 1-loop computation of the wave function corrections of 
A , Ai and the gauge coupling g 2 . We proceed as explained in [|J, changing pf — > Pf—ifi 
in the fermion propagators, and including the 2-loop contributions into the effective 



potential from [43|, jig] . Let 

1 



V(u) = (4tt) 



{(pf + inTujy + p 2 ) 2 (p} + p 2 ) 2 



j Pf l y \p f t *vi ± lu) t p ) ypf 
dp/l-e p 1 1 



00 p V 1 + eP eP +nuj + 1 e~P+™ + 1 



(5.4) 
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At small u, V{u) « -(7C(3)/2)o; 2 . We then obtain 



2 2 
03= 03 



ft=0 



iV 



6 + N c — Nf 
N f 



r- r \ I 1 ©(-^ - ^' 2 

' x U=ol i 6 + iV c -JV> W 



a; + C(x 



VN C =2 : 



8 _ i y / )( 4 + Ar / ) / 3 iV> , fi 
U4n 2 x \ 4+IV,WT 



4 + JV> v ?rT 



192 - 2N f - 7Nj - 2Nf / 3^(11 + 2N f )(A - N f ) , fi ^ 
Wtt2(K - NA V + 192 - 2N f - 7N 2 - 2Nj Wt) 



VN C =3 



96vr 2 (8 - N f 
(9-iV»(6 + JV» 



1447r 2 x 



1 + 



3iV / ( V \ 2 



+ 



6 + N f \nT' 
1 + 



486 - 33N/ - LUV) - 2AT) / 3Nf{7 + N f )(9-2N f ) , \l x2 
96vr 2 (9 - A/) V 1 + 486 - 33 N f -llN 2 -2N 3 f WT' 



u N f / 3 7 + Nt n, 
z Nc=3 = £^(1 + -- ^-x + 0{x\ 



T 3vr 2 V 29- N f 
In the case we actually study, N c = 3, Nf = 2, these reduce to 



(5.5) 
(5.6) 
(5.7) 

0(x), 

(5.8) 

+ 0(x), 
(5.9) 



x 



xi,._^l- 2 -V(^ ]X 



\fj,=0 



7 15 \/ 3/ji_N2\ 

18vr 2 x + 28tt 2 J V + A^tiT> )' 
fi 2 



, 27 

, 1 H x 

ttT 3tt V 14 



(5.10) 
(5.11) 
(5.12) 
(5.13) 



where gf , 



x 



jU=0 



are from Eq. 



and corrections are of relative order 0(x 2 ). Let us 



recall that x ~ 0.1 at T ~ 2T C . 

Due to the small value of x, we will in the analysis which follows for simplicity 
ignore even the 0{x) corrections in g%,x,z in Eqs. (|5.10|) , ( |5.11|) , ( [5.131 ), and use the 
leading order expressions displayed already in Eq. ( |5.3| ) . All of the sub-dominant effects 
neglected work in the same direction and would strengthen the changes we see in the 
mass spectrum. For z this is obvious because, according to Eq. ( |5.13| ), for a given 
fi/T the physical value would be larger than what we have used. For x, Eq. ([5.11 ) 



and the negative value of T> imply that the fact that we have kept x unchanged means 
that we have effectively moved slightly up in T/A^g. In g 2 a major part of this effect 
cancels since the correction is the same in Eqs. ( |5.10| ), ( |5.11| ). Therefore, M/T remains 
essentially the same for given M/g 2 , even though the temperature is slightly higher. 
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5.2 A complex action and reweighting 

The generic problem of SU(3) simulations at non-vanishing chemical potential now 
appears as follows. For /i 7^ 0, there is the imaginary term izS z in the action, and 
the conventional importance sampling will not workf]. A way to circumvent this is to 
do the importance sampling with the action at z = 0, while exp(— izS z ) is included 
in the operators. We will refer to this procedure as "reweighting". Reweighting is 
also applied in the "Glasgow method" suggested for ji 7^ simulations in 4d (for a 
review, see f46|); however, in our approach we do not need to carry out any expansion 
related to z. Denoting by brackets the expectation value with respect to the action in 
Eq. (|3.1| ), operators even + and odd 0„ in R therefore appear as 



(0 + exp(-izS z )) = (0 + cos(zS z )), (CL exp(-izS z )) = -i(O-W0.(zS z )). (5.14) 



Since Eq. (|5.14 ) involves trigonometric functions, the problem manifests itself if con- 



figurations with S z ^> 1 occur frequently. In this case there are cancellations in the 
functional integral, and the Monte Carlo method will soon lose its accuracy. 

Configurations with S z 1 need not always be typical, however. Indeed, let us 
consider the distribution of S z in Eq. ( |5.2| ) for z = 0. The centre of the distribution is 
at zero. The fact that there is Debye screening, y > 0, guarantees that the distribution 
is approximately Gaussian. Its width is easily computed, and is given by the 2-loop 
scalar sunset graph on the lattice, which has been evaluated in pE |26| : 



Z /6L\3/2 



A / m \ Z /ULj 



M 1/2 
• r ) ( In + 0.58849 * 



(5.15) 



We note that the width grows with the lattice size L (or, more precisely, with the 
physical lattice size L//3)F1. Nevertheless, if z is small enough, the infinite volume and 
continuum scaling limits of the correlation functions may be extracted before Eq. ( [5. 15] ) 
grows to values 3> 1. In these cases we can in practice carry out all the measurements 
before the oscillations set in. The feasibility of this procedure can be checked by 
monitoring the numerical distribution of zS z , the argument of the reweighting factors. 
Accurate calculations should be possible whenever the bulk of the distribution is con- 
tained in the interval (— tt, tt). Once the distribution gets significantly broader, the 
breakdown should show up in the correlation functions simply as a noisy signal. 

Even assuming that the oscillations are under control, reweighting according to 
Eq. (|5.14j) of course requires some care. The importance sampling of the Monte Carlo 



is dominated by the minimum of the action at z — 0, which is unmodified by including 
S z in the operators. If the z = minimum is widely separated from the true 2^0 min- 
imum in configuration space, there is a danger of biasing the Monte Carlo towards the 



4 It can be seen directly in 4d that the sign problem is strongly correlated with the imaginary part 
of the temporal Polyakov loop This corresponds precisely to zS z in the 3d language. 
5 For a fixed physical lattice size, the width grows with f3, but only logarithmically. 



18 



wrong configurations. In order to check this, we have tested the reweighting method 
by considering z imaginary, z —> In this case the action is real and the scalar 

update can be modified by a Metropolis step to include S z , and hence update with the 
exact action. The results can then be compared with those obtained from Eq. ( |5.14j ), 
now with 

cos(zS z ) — > cosh(£S z ), — ism(zS z ) — > — sinh^S^). (5.16) 

Simultaneously, we change the sign of /i 2 in y in Eq. ( |5.3|) . 

Finally, as a small technical detail let us note that in practice we implement S z in 
Eq. ( |5.2| ) on the lattice by including an 0(a) improved []27| version of TtAq, as in 
Eq. (2.10) of P|. 



5.3 Determining masses with a complex action 

Assuming that our reweighting procedure works, we then have to extract masses in the 
changed symmetry situation. After the inclusion of /x ^ 0, the action of the effective 
theory is no longer real, and no longer invariant in R, C . It is however still invariant 
under RC and P. In the scalar channel this means that the operators and 0t~ in 
Eq. fl2.13|) can couple to each other, and similarly and CT + . These two channels 



are still distinguished, however, by the parity P. 

We measure the correlation matrix between the two channels which used to be de- 
coupled for fj, = 0. From Eq. (|5.14 ) we know that the result is of the block form 



where A(t), B(t) are symmetric matrices (possibly of different sizes), representing the 
correlations within the two channels as in Eq. ( |3.6| ). We have factored out z in the 
off-diagonal blocks, to make it clear that the result is odd in z and vanishes for z — > 0. 
In order to facilitate the solution of the eigenvalue problem in Eq. (|3.7| ), we may note 
that we can equivalently consider the eigenvalue problem for the real correlation matrix 



C(t) = u,X> , (5-18) 





obtained from C(t) by a similarity transformation Z = diag(l,z). Finally, for imagi- 
nary chemical potentials z = —i( we have a correlation matrix of the usual symmetric 
form as in Eq. (|3.6|), 

In addition, let us note that physical observables such as masses (i.e., inverse corre- 
lation lengths) must be even in z, since a change of its sign can be compensated for by 
the field redefinition A Q — > —A Q in the action in Eq. ( |5.3|) . Since the mass spectrum 
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is well defined and there are no massless modes at z = 0, we moreover expect that 
the system is analytic in z for small \z\, the expressions of the masses starting with 
z 2 . In particular, masses must remain real for small enough z (even though for the 
eigenvalues of a fully general matrix of the form in Eqs. fl5.17| ), ( |5.18| ) this is not always 
the case). 

We can use these observations to discuss the form of the mass eigenstates. For the 
cases in Eqs. ( |5.17|) , ( 15.18 ), we do not in general have orthonormal eigenvectors in the 



usual sense. In order for the procedure to make physical sense, we however expect 
(and observe) that the eigenvectors are independent of t within statistical accuracy, 
with eigenvalues of the form ~ exp(— aMt), as in Eq. ( |3.7|) . 

Consider now first Eq. ( |5.19| ). Since the eigenvalues are even in (, the eigenvectors 



should be of the forms 

v (1) = ( c V v A b ) , vU(^ = l; v(2) =("vr)' ^ + v ^ = 1 ' ( 5 - 2 °) 

where normalisation conditions have also been shown. With analytic continuation, we 
can then directly apply this to the case of Eq. ( |5.17|) , 



vi-^ 1; v«^-;;-j, -,M + vi = l, (5.2!) 
as well as to the case of Eq. ( |5.18j ), using the similarity transformation: 



v« = { ^ b j , vi - z\l = 1; = { j , -z\\ + v| = 1. (5.22) 

We observe that the "metric" has changed, but otherwise the situation is quite analo- 
gous to the usual one in Eq. (|5.20|) . 



Fortunately, in practice even less needs to be changed in the usual numerical analysis 
based on a matrix of the form in Eqs. (|3.6|) , ( |5.19| ), if we use Eq. ( |5.18|) . There 
are algorithms finding the eigenvalues and eigenvectors for a general non-symmetric 
real matrix. The only change is in the normalisation according to Eq. ( [5.22| ), which 
amounts to fixing an overall coefficient for each eigenvector. However, as we only 
consider the correlations between the eigenvectors, cf. Eq. ( |3.8| ), and do not need to 
verify orthogonality explicitly, we can equally well employ a "wrong" scalar product 
and normalisation based on the old type of metric in Eq. ( |5.2(J| ) (with £ — > z). Then 
everything goes precisely as before. 

Finally, let us illustrate the general expectations for the mass pattern with a simple 
2x2 matrix as an analogue of Eq. (|5.17p . Consider two real scalar fields 0, x, with a 



mass term 



V = ^m 2 4> 2 + X -M 2 X 2 + ie^x- (5.23) 



20 



fi/T 


y 


z 


L 2 ■ 


T 


0.5 i 


0.47382 


0.0338 i 


20 3 




1 n ; 

l.U 1 




U.UO 1 O 1 


on 3 




1.5 i 


0.40042 


0.1013 i 


30 3 




0.5 


0.49218 


0.0338 


30 3 




1.0 


0.51970 


0.0675 


30 3 




1.5 


0.56558 


0.1013 


30 3 




2.0 


0.62981 


0.1351 


30 3 , 


40 3 


4.0 


1.07026 


0.2702 


10 3 , 


14 3 ,18 3 ,30 3 



Table 4: The lattice parameters and sizes used for calculations with // ^ 0. All are for 
T = 2A M g, N f = 2, x = 0.0919, (3 = 21. The parameter z is from Eq. fi5l\) . 



For M ^> m and e < (M 2 — m 2 )/2, the mass eigenvalues are real, m 2 + (e/M) 2 , 
M 2 - (e/M) 2 . Thus the lightest mass goes up, which is what we expect for a real 
chemical potential. In fact there is another effect contributing in the same direction, 
since according to Eq. ( |5.3| ) the parameter y grows with \xjT . In the case of a complex 
e = —17], on the contrary, the lightest mass becomes even lighter, the heaviest mass 
even heavier: m 2 - (r//M) 2 , M 2 + (r//M) 2 . A gain, the change of y contributes in the 
same direction for the lightest mass. 



6 Numerical results for /i ^ 

We now proceed to present the results of our first numerical explorations of fi ^ 0. We 
restrict our attention to the case T = 2T C , Nf = 2, corresponding to x = 0.0919. As 
we have explicitly checked for good scaling behaviour at fi = 0, we work exclusively 
with (3 = 21 in this section. The parameters and lattices considered are summarised 
in Table 0, and the detailed results are displayed in Tables ITT], O in the Appendix. 



6.1 Consistency checks and the general pattern 

We begin by assessing the distribution of the argument of the reweighting factors in 
Eq. ( |5.14p . We find that any typical histogram can be easily fitted to a Gaussian, from 
which we determine its width. The plots in Fig. |5| show the corresponding distributions 
as a function of fi/T and the lattice volume, with behaviour as expected from Eq. (|5.15|) . 
The volume chosen for the fi/T series is large enough to be free from finite size effects in 
the n = simulations. As the histograms show, for these volumes the distribution of the 
argument of the reweighting factor is entirely contained within one period for fi/T < 
2.0. On the other hand, at fi/T = 4.0 the tails of the distribution are significantly 
spreading. Nevertheless, the bulk of the distribution is still within the region without 
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Figure 5: Distribution of zS z [A ] for the reweighting procedure (Eq. ( \5.14\ ))- 



sign flips. For the same number of measurements at /i/T = 2.0, 4.0, however, the errors 
on the correlation functions are indeed by a factor of three larger in the latter case, 
confirming the onset of the sign problem. Going to still larger fi/T will thus lead to 
increasingly noisy signals, until a mass determination becomes impossible. 

In the next step, we test the reliability of the reweighting procedure by comparing 
with the correct importance sampling in the case of an imaginary chemical potential. 
This comparison can only be carried out for a limited range of /a/T. In the full action, 
the imaginary chemical potential term weakens the metastability of the symmet- 
ric phase, and tunnelling during typical simulations to the unphysical phase becomes 
possible. This can be countered by increasing the lattice volume, but in so doing we 
increase the width of the distribution of z^LAo]. We find however that there exists a 
window of opportunity between these two limits where we can usefully compare the 
two algorithms. In the left panels of Fig. [| the ground states of the + and 0~ chan- 
nels are shown as a function of imaginary /i/T. First, we note the correct qualitative 
behaviour of a mass decrease with \(J,/T\, as described after Eq. ( |5.23| ). Second, we 
observe complete agreement within the small statistical errors between the reweighted 
observables and those measured with the full action. 

With these tests passed, we can thus display the same states for the case of real 
chemical potential in the right panels of Fig. || Again in agreement with the qualitative 
behaviour expected from the discussion after Eq. (|5.23|) , the lowest mass values are 
growing with increasing /i/T, since the constituent mass parameter y 1 / 2 gets larger, 
see Eq. (|5.3| ). Thus the corresponding correlation lengths decrease. We note again, 
however, that purely gluonic states (or, in the 4d language, purely magnetic states 
~ Tr F 2 2 ) are rather insensitive to /i/T, in analogy with their behaviour under variations 
of y discussed earlier. This means that at large enough /i/T > 4.0, they may become 
the lightest states in the system, as suggested by the top right panel of Fig. |6|. 

At fi/T = 2.0 we have performed a finite volume check, comparing lattices of L = 30 
and 40 at (3 — 21. We find the low lying masses to be consistent within statistical 



errors, see Table 12. This is encouraging, implying that the reweighting does not 
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Figure 6: TTie dependence of the lowest thermal masses on fi/T for imaginary (left) 
and real (right) chemical potential. For clarity of presentation, we show the lightest 
state only for each quantum number assignment. The quantum number assignments 
refer to the eigenstates at fi/T = 0. Filled symbols correspond to glueballs. "Exact" 
and "reweighted" simulations (left) are explained in Sec. \5.3j . 



magnify finite volume effects unduly, and that it is indeed possible to extrapolate 
screening masses to the infinite volume in a finite chemical potential context, before 
the onset of oscillations. 



6.2 The Debye screening length at finite baryon density 

Let us now consider what happens to the Debye screening length as defined in Sec. [4 . 1 . 2| . 
This definition is based on the i?-symmetry of our effective theory pi . For /i ^ the 
extra term in the action spoils the symmetry, resulting in a merging of the various 
.RC-channels, while only P remains as a good quantum number. The question then 
arises, how can the Debye mass be defined in such a situation [24|? 

We now note the following. For \i = 0, the lightest R = — 1 state was 0z + ~ Tr AqF u . 
This is, however, also the lightest state in the channel P = — 1. Thus, even after /i ^ 0, 
the same operator continues to determine the long distance decay of a number of 
R = — 1 operators. The values can be seen in Fig. |6], bottom right. If we define this 
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largest P = — 1 correlation length to correspond to Debye screening, then fi ^ does 
not change the situation in an essential way, and the corresponding correlation length 
decreases with \i. An example of a 4d gauge invariant operator with such a behaviour 



isTrF 03 F 12 @. 

On the other hand, the operator ~ Tr A\ does couple to 0+ + ~ Tr after the 
inclusion of 7^ 0. Therefore, its correlation length strictly speaking increases by a 
factor of almost two as soon as fi 7^ 0, as can be seen from the top right panel in Fig. || 
There we have labelled the states with the quantum numbers of the fi/T = limit; 
however, for /i/T > 0, any operator which originally had just the quantum numbers 
0l~ such as Tr Aq, couples now also to 0^ + . The eigenstate shown in Fig. |^ which 
does have a larger mass, is a particular linear combination of the two types of original 
states. Nevertheless, the overlap between the two operators ~ Tr^Q, ~ TrAjj is very 
small for small fi/T, so that even for ~ TrAjj there is an intermediate distance range 



where correlations should decay as shown with 1 in Fig. 0, as discussed in [24 . Thus 



no abrupt change is observed. An example of a gauge invariant 4d operator with such 
a behaviour is the imaginary part of the 4d temporal Polyakov loop. 



7 Conclusions 

In this paper, we have carried out mass measurements in the different quantum number 
channels within the 3d SU(3)+adjoint Higgs theory. We have stayed on the dimensional 
reduction curve in the symmetry restored phase of the 3d theory, and the results are 
thus expected to correspond to spatial correlation lengths in the deconfmed quark- 
gluon plasma phase of QCD. This interpretation appears to be relatively accurate at 
least down to T ~ 2T C . We first considered a vanishing chemical potential [i = 0, and 
then also extended the measurements to a finite fi/T< 4.0. 

For n — 0, Nf = 0, we believe that the asymptotic correlation lengths related to 
the real and imaginary parts of the 4d temporal Polyakov loop, as well as to other 
gauge invariant bosonic operators, are now relatively well understood. For instance, 
at T ~ 2T C « 2Aj^g, the real part of the Polyakov loop (with quantum numbers 
according to our conventions discussed in Sec. |2TT| ) decays exponentially at ~ (3T) _1 , 
while the imaginary part (0l~) decays at ~ (5T) _1 , as can be observed from Fig. III. The 
Debye screening length according to the definition of [0] turns out to be determined 
by operators of the type CT + ~ TrF 3-F 12 , and is slightly longer than the screening 
length of the imaginary part of the Polyakov loop, ~ (4.6T) -1 . 

We have also explicitly studied the effects of Nf dynamical fermions on the longest 
correlation lengths. For fixed T/A^jg, the correlation lengths decrease as Nf is in- 
creased. For instance, in the phenomenologically interesting case of Nf = 2, 3 at 
T ~ 2A^jg, the real part of the Polyakov loop decays exponentially at ~ (4T)" 1 , while 
the imaginary part at ~ (6T) _1 . As a theoretical point, it is interesting to note that 
in units of gf , the screening lengths scale quite well with A^ c . 
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We have then extended the measurements to /x 7^ 0. We have demonstrated that the 
use of dimensional reduction allows one to carry out simulations corresponding to phe- 
nomeno logically interesting values of quark masses, temperatures, and /x/T (although 
not at the point of the phase transition). Simulations are possible because the only 
term in the super-renormalizable 3d action suffering from the sign problem, ~ iTi Aq, 
comes with a very small numerical coefficient. 

In general, we observe that as /x/T is switched on at a fixed temperature, the screen- 
ing lengths decrease. This can be interpreted so that the lightest excitations, involving 
Aq, are "less critical" and correspondingly heavier. Thus, staying at a fixed temperature 
but increasing /x, we are apparently moving further away from the critical line, which 
hence has to bend down. This is in complete accordance with what we qualitatively 
know about the phase diagram in the (/x, T)-plane. 

However, there is one state whose mass does not increase, the 3d glueball ~ Tr F^ 2 - 
Thus at large enough /x/T, similarly to the case of large enough T/Aj^g at /x = 0, it 
becomes the lightest excitation in the system. We find that the crossover corresponds 
roughly to y ~ 1.1 in terms of the dimensionless parameter defined in Eqs. (|2.2j) , (|5.3|) . 
For Nf = 2, this corresponds roughly to T/Aj^g ~ 10 2 at /x/T ~ 0.0, and /x/T ~ 4.0 
atT/A^-2. 

Finally we have addressed the behaviour of the imaginary part of the 4d temporal 
Polyakov loop. As the chemical potential is switched on, its asymptotic correlation 
length suddenly increases contrary to the other observables discussed, and agrees for 
all /x ^ with that of the real part of the 4d temporal Polyakov loop. 

The task remains to establish a direct connection between the screening masses 
discussed here and observables relevant for the phenomenology of heavy ion collision 
experiments. 
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Note added 

After the submission of our paper, a letter appeared by Gavai and Gupta |fT7| , who 
measure the spatial meson (pion) correlation length at T ~ (1.5...2.9)T C for Nf = 4 
light dynamical fermions. They observe that at T ~ 2T C the meson "mass" measured 
from ipif) has become lighter than the lightest bosonic mass (represented by ~ Tr Aq 
in the 3d language), so that the bosonic gauge degrees of freedom are no longer the 
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lightest dynamical degrees of freedom, and thus the standard dimensionally reduced 
effective theory should not be reliable any more. This observation is in accordance with 
our comments in Sec. |4.1.1| . However, Gavai and Gupta make the further suggestion 
that even in the temperature range T ~ (2...10)T C where the meson mass is no longer 
the lightest one, but is still below its asymptotic value (~ 2ttT in the case of the 
chiral continuum limit), dimensional reduction may not be reliable, in contrast to the 
situation in the purely bosonic case Nf = 0. Let us comment that there is no hard 
evidence for the latter suggestion. In fact, it seems quite possible to us that the decrease 
from ~ 2nT could be partly accounted for precisely by the confining bosonic 3d gauge 
field dynamics related to Aq, A^. Furthermore, we see no evidence for why the bosonic 
correlation lengths we have measured here would not be satisfactorily reproduced by 
the 3d theory in this temperature range. Thus we believe that our bosonic mass 
measurement for Nf > (both at fi/T = and /i/T > 0) do reproduce at least the 
qualitative features of the physical 4d theory at all temperatures T > 2T C . 

As already discussed by Gavai and Gupta, one might also expect that chiral effects 
responsible for the behaviour seen by them are at Nf = 4 somewhat larger than in 
the realistic case Nf = 2. ..3. For instance, Tr is lighter for smaller Nf as we have 
discussed in this paper, while ipip still reaches the same asymptotic value, so that 
the level crossing invalidating dimensional reduction should take place at a smaller 
temperature. 
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Appendix A Tables 

For completeness, we collect in this Appendix the numerical results for the correlation 
length measurements we have carried out. 
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Table 5: Mass estimates and dominant operator contributions in the 0+ + and 2^ + 
channels for Nj = 0. The dominant operator types contributing are denoted with 
(without) parentheses if {<&\<f) k ) < (>)0.5. 
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Table 6: Mass estimates and dominant operator contributions in the 0+ + and 2+ + 
channels for various Nf > on L = 30, 40 for (3 = 21, 28 respectively. The dominant 
operator types contributing are shown in the superscripts. 
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Table 7: Masses M/gl for channels other than 0\ + ,2~\_ + at Nf = 0. The stars in the 
superscripts denote excited states. 
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2.83 (7) 




3.71 (21) 


3.94 (19) 


3.70 (18) 


3.58 (28) 


3.84 (25) 


3.64 (14) 


o:+ 


1.932 (32) 


1.977 (43) 


1.960 (43) 


2.079 (35) 


2.103 (67) 


1.932 (43) 


o:+* 


3.46 (11) 


3.52 (21) 


3.50 (21) 


3.56 (11) 


3.66 (16) 


3.56 (14) 


01" 


2.19 (7) 


2.27 (11) 


2.391 (35) 


2.604 (70) 


2.677 (70) 






3.31 (24) 


3.29 (19) 


3.31 (15) 


3.53 (21) 


3.60 (11) 


3.29 (10) 


1+ 


3.59 (25) 


3.92 (20) 


4.06 (16) 


4.11 (32) 


4.08 (25) 


3.73 (15) 


1_ 


3.10 (18) 


3.17 (25) 


2.91 (16) 


3.28 (28) 


3.21 (18) 


3.08 (12) 



Table 8: Masses Mjg\ for channels other than 0| + ,2| + ; for various Nf > 0, on 
L = 30, 40 with (3 = 21, 28, respectively. The stars in the superscripts denote excited 
states. 





T = 2T C 


T ~ 10 n T c 


(3 = 21 


(3 = 28 


(3 = 28 


L = 38 


L = 28 


L = 40 


L = 40 


aMp(L) 


1.042 (40) 
0.167 (4) 
0.585 (15) 


0.749 (20) 
0.166 (3) 
0.580 (11) 


0.585 (17) 
0.122 (2) 
0.569 (10) 


0.625 (20) 
0.126 (2) 
0.588 (10) 



Table 9: Spatial Polyakov loop masses and string tensions for Nf = 0. 





T = 1.5A m 


T = 2.0A^ 


N f = 2 


N f = 3 


N f = 2 


N f = 3 


N f = 4 


N f = 2 


(3 = 21 


(3 = 21 


(3 = 28 


aMp(L) 


0.760(21) 

0.161(3) 

0.563(11) 


0.818(15) 

0.167(2) 

0.584(8) 


0.798(19) 

0.165(2) 

0.577(8) 


0.847(35) 

0.170(4) 

0.595(15) 


0.820(19) 

0.167(2) 

0.584(8) 


0.588(11) 

0.123(2) 

0.574(10) 



Table 10: Spatial Polyakov loop masses and string tensions corresponding to 3d pa- 
rameters with Nf > 0, on L = 30, 40 for (3 = 21, 28 respectively. 
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/i/T = 0.5i 


li/T-- 


- l.Oi 


ii/T-- 


= 1.5i 




L = 20 


L = 20 


L = 20 


L = 30 


L = 30 




rewtd. 


rewtd. 


exact 


rewtd. 


exact 


0+ 




M/gl 


M/gl 


M/gl 


M/gl 




1.34 (3) 


1.20 (2) 


1.18 (2) 


0.89 (3) 


0.85 (3) 




1.94 (6) 


1.89 (6) 


1.73 (11) 


2.54 (14) 


2.96 (21) 




2.21 (9) 


2.22 (7) 


2.21 (11) 


1.81 (6) 


1.81 (6) 


$f 


2.54 (18) 


2.53 (9) 


2.50 (13) 


2.45 (11) 


2.45 (18) 


o- 


M/gl 


M/gl 


M/gl 


M/gl 


M/gl 




2.04 (3) 


1.87 (4) 


1.86 (3) 


1.76 (3) 


1.77 (7) 




2.91 (15) 


2.88 (9) 


2.81 (13) 


2.46 (7) 


2.42 (9) 




3.15 (25) 


3.32 (18) 


3.45 (20) 


3.29 (21) 


3.23 (25) 




3.80 (35) 


3.75 (25) 


3.63 (28) 


3.50 (27) 


3.70 (28) 



Table 11: Mass estimates and dominant operator contributions in the + and 0~ chan- 
nels for imaginary /i ^ from reweighted and exact actions. The dominant operator 
types contributing are shown in the superscripts. Some states on the lattice size L = 30 
(used for /i/T = 1.5i) have been reordered for clarity of presentation. 
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ii IT - 


-05 


ii IT - 


- 1 


ii IT - 


- 1 Fi 




L = 


30 


L = 


30 


L = 


30 


0+ 


M/gl 


Ops. 


M/g'i 


Ops. 


M/gi 


Ops. 


$1 


1.37 (4) 


R 2 , L 


1.48 (2) 


R 2 , L 


1.59 (2) 


R 2 , L 




2.29 (8) 


C, (R 3 ) 


2.61 (3) 


C, (R 3 ) 


2.50 (5) 


C 


$ 3 


2.51 (9) 


Rs, (C) 


2.61 (4) 


Rs, (C) 


2.78 (5) 


R 3 


$4 


2.64 (17) 


P 


2.84 (8) 


P 


2.82 (7) 


P 


0" 


M/g'i 


Ops. 


M/g'i 


Ops. 


M/gi 


Ops. 


$1 


2.01 (2) 


B l 


2.02 (2) 


B 1 


2.14 (3) 


B l 


$ 2 


2.84 (12) 


B 2 


3.12 (6) 


B 2 


3.18 (7) 


B 2 


$ 3 


3.49 (18) 


B l 


3.56 (7) 


B 1 


3.67 (11) 


B l 


$4 


3.74 (21) 


C 


3.68 (11) 


C 


3.73 (12) 


C 









= 2.0 




Vi/T = 


: 4.0 




L = 


30 


L = 


40 


L = 


30 


0+ 


M/g'i 


Ops. 


M/g'i 


Ops. 


M/gi 


Ops. 


$i 


1.73 (3) 


R 2 , L 


1.80 (11) 


R 2 , L 


2.58 (16) 


R 2 , L 




2.45 (8) 


C 


2.43 (21) 


C 


2.69 (18) 


C 


$ 3 


3.01 (14) 


Rs 


2.71 (21) 


Rs 






$4 


3.02 (14) 


P 










0" 


M/gl 


Ops. 


M/g'i 


Ops. 


M/g'i 


Ops. 


$1 


2.25 (5) 


#i 


2.19 (7) 


B l 


2.43 (20) 


#i 


$2 


3.45 (8) 


B 2 


3.37 (25) 


B 2 


4.17 (22) 


B 2 ,C 


$ 3 


3.57 (11) 


C 


3.78 (35) 


C 






$4 


3.70 (16) 


B 1 


3.85 (36) 


B 1 







Table 12: Mass estimates and dominant operator contributions in the + and 0~ 
channels for real /i ^ 0. The dominant operator types contributing are denoted with 
(without) parentheses if \ |< (>)0.5. 
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